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Abstract 

The Hamiltonian formulation of superfiuids based on noncanonical Poisson brackets is studied 
in detail. The assumption that the momentum density is proportional to the flow of the conserved 
energy is shown to lead to the covariant relativistic theory previously suggested by Khalatnikov, 
Lebedev and Carter, and some potentials in this theory are given explicitly. We discuss hydro- 
dynamic fluctuations in the presence of dissipative effects and we derive the corresponding set of 
hydrodynamic correlation functions. Kubo relations for the transport coefficients are obtained. 
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I. INTRODUCTION 



Recently there has been considerable interest in the possible occurrence of superfluidity 
in relativistic systems such as neutron stars or quark matter at very high densities [1]. A 
common distinctive feature of these systems is the spontaneous breaking of (at least) a U{1) 
symmetry associated with particle number. 

At nonzero temperature, the dynamic description of a given system at large distance and 
time scales is based on hydrodynamic equations for the hydrodynamic modes. Generically 
these consist of quantities whose long-wavelenght fluctuations have a large lifetime which 
becomes inflnite as k — 0. The hydrodynamic variables include densities of conserved 
quantities and Goldstone modes of broken symmetries. 



A relativistic extension of superfluid hydrodynamics of "^He [2|, [l2| was previously sug- 
gested by Khalatnikov, Lebedev and Carter in Refs. jl, 0,[5j. More recently, Son outlined the 
construction of the nondissipative hydrodynamics of relativistic systems with broken sym- 
metries using the Poisson bracket method in the case of broken U{1) symmetry ^ and in the 
case of nuclear matter [71], where the chiral symmetry SU (A^/)l x SU (A^/)_r is approximately 
broken down to SU{Nf)L^fi. 

In this article we will pursue the scarce research on this subject. In particular, in Sec. [TTl 
we derive in detail the equations of motion of superfluids using the Poisson bracket method, 
and we show that the assumption that the momentum density is proportional to the flow 
of the conserved energy leads to relativistic hydrodynamics of Khalatnikov and Lebedev. In 
this framework, with the alternative assumption that the momentum density is proportional 
to the flow of the U{1) conserved charge one obtains the two-fluid model of non-relativistic 
superfluidity. In Sec. IIIII we quickly review the memory function formalism to study the 
hydrodynamic fluctuations and examine some sum rules, following closely the treatment 
given in Ref. [l4|. We use this information, together the information about the forces and 
the memory matrix, to deduce the correlation funcions, the Kubo relations and the linearized 
constitutive relations in Sec. [IVl As well, some Ward identities are checked. 

Although perhaps none of the statements made below could go beyond what is known 
about the equations of motion of relativistic superfluids, it can be of some interest to present 
the discussion of the main results from the point of view of the theory of hydrodynamic 
fluctuations and correlation functions, in order to show how this general framework organizes 
the physical description in the low-energy limit. 



II. POISSON BRACKET AND HYDRODYNAMIC EQUATIONS 

We flrst consider the low energy, nondissipative dynamics of a system without long-range 
interactions with broken U{1) symmetry at a temperature well below criticality. In such a 
system, the hydrodynamic variables are flve conserved densities, namely, the entropy per 
unit volume s, the density n of the U{1) charge and the momentum densities g^, plus a 
Goldstone mode (p. The dynamics is governed by the Hamiltonian functional H[s, n, g\ djip] 
whose speciflc form can be computed or guessed from the underlying microscopic model. 
Note that the invariance of H under U{1) transformations, 6ip = a, prevents the dependence 
upon (p. The energy functional is an extensive quantity 

H = I d^xe, (2.1) 
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where e{s,n, g^,dj!f) is the energy density from which one can compute thermodynamic 
equihbrium properties. As the following discussion concerns the hydrodynamic behaviour 
for large wavelength and low frequency, it will be enough to consider the dependence of e 
up to first derivatives of ip. In order to obtain all the equations of motion, we must specify 
a Poisson bracket structure [F, G]pb between functionals of the hydrodynamic variables. 
Then, the time derivative of a functional V of {s, n, g^, (p} can be derived from 

dV 

^ = [V,^]PB- (2.2) 

The noncanonical Poisson bracket is taken to be 

'6F d 6G 5G d 5F 



+ n 
+ s 



Sg^ dx^ Sg^ Sg^ dx^ 5g^ 
6£_d_6G _ 6G d 6F' 
5g^ dx^ 6n 6g^ dx^ 5n 
6F d 6G 6G d 6F 



6g^ dx^ 6 s 6g^ dx^ Ss 
/6F6G 6G6F\ f5F5G 5G 5F\ 



dx^ \Sg^ 5(p 5g^ 5(p J \5n Sep 5n S^p J 



(2.3) 



The form of the bracket in the first three terms follows from the conservation of the linear 
momentum, the f/(l)-charge and the entropy The remainder terms refiect the transfor- 
mation property of the (^-field under infinitesimal spatial translations and the fact that the 
charge density and the Goldstone mode are canonically conjugated 0] 

Kx),^(y)U = 5(x-y). (2.4) 

With this normalization the Goldstone mode ip is dimensionless. 

In terms of the quantities conjugate to the hydrodynamic variables, the temperature T, 
the chemical potential /i, the fiuid velocity v and the vectorial quantity A = de/dCVip), the 
energy density e satisfies the thermodynamic relation 

de = Tds + fidn + v ■ dg + X- d{Vip). (2.5) 

A second thermodynamic identity to be used below is the Gibbs-Duhem relation for the 
pressure, defined by 

p = -£ + Ts + /in + V ■ g. (2.6) 

The non occurrence of A ■ V<y9 in this expression can be understood by noting that the 
hydrodynamic Goldstone mode, which is not a conserved density, cannot contribute to the 
total (integrated) entropy S since it corresponds to a single coherent mode [l^. Thus, the 
entropy is a homogeneous function of all extensive variables: the volume V and the conserved 
quantities, energy, linear momentum and U{1) charge. Accordingly, 

^ dS^^ dS ^ ^ dS : dS ^ 

i 

which, after multiplication by V^'^, yields the Gibbs-Duhem relation. Notice, however, that 
the thermodynamic identity for the pressure has the form 

dp = sdT + nd^ + g-dv- X-d{V(p). (2.8) 



3 



From these formulas and the Poisson bracket it is easy to derive the equations of fluid 
dynamics without dissipation. Let us list them. The conservation laws adopt the form 

9ts(x,t) + V-Gsv) = 0, (2.9) 
(9tn(x,t) + V- (nv + A) =0, (2.10) 

dtg\x,t) + dkt'' = 0, (2.11) 

where the reactive part of the stress tensor is given by 

t''' =pS''' + v''g' + X'^diifi, (2.12) 
and the equation of motion for (f is 

dt(f{x,t) = -fi-v -Vip. (2.13) 

As a consequence of these equations of motion and the Gibbs-Duhem relations (12. 6p 
and (12.81) . we find the flow of the conserved energy 

dt8{^,t) + dal = 0, (2.14) 
ji = {f, + ^.V^)y + {e + py. (2.15) 

At this point, we can separate the momentum density g and A into two pieces proportional 
to V and Vyj by-*^ 

g = aiv + q;2Vv9, (2.16) 
A = asv + ^Vy?, (2.17) 

where the a^,^ are considered as functions of {s,n,v\dj(t>). The quantity ^ has dimen- 
sion mass^ and, as we shall see below, characterizes the long wave limit of the equal-time 
correlator of (p. The symmetry of t''^ requires that = —ai- 

Further progress can be made if we enforce Lorentz invariance of the hydrodynamic 
equations. This can be made by assuming that the momentum density is the same as the 
energy flow, j* = . Consequently, from Eqs. (12.151) and (I2.16p . we find 

ai = e + p- e(/i + v- Vipf , (2.18) 
"2 = e(/i + v- V(^), (2.19) 

which lead to the constitutive relations 

f J = + [5 + p _ ^ (^ + V ■ V(^)^] t;V + idi^pdj^p, (2.20) 
i = [£ + p-e(/i + v- V<^)>* + e(/i + v- V^)9,<^, (2.21) 
A' = -e (/i + V ■ V^) v' + id,^. (2.22) 

It is important to emphasize that the energy flow is itself conserved because of the assumption 
of relativistic invariance. This fact will have implications on the form in which the dissipative 
terms must be included in the full constitutive relations. 



^ Because of rotational invariance, this decomposition is completely general and does not amount to a loss 
of generality. 
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With the introduction of the four-vectors 



^^ip={-^Ji-^,.Vip,Vip), (2.23) 

= (T + v w,-w), (2.24) 

J'^ = (n,nv + A), (2.25) 

s^ = (s,sv), (2.26) 

where w = s~^(g — n'Vip), it turns out that the differential Gibbs-Duhem identity can be 
written in a covariant form as dp = — J'^d{d^(p) + s'^dWfj,, and the energy-momentum densities 
and the stress tensor, e = T^^,jl = T^^,f^ = T*-', are obtainable from the symmetric four- 
tensor 

Tl^u ^ _Jt^QU^ ^ g^^u _ ^f.u_ (2.27) 

Correspondingly, the conservation equations can be written as 

d.T^^" = s'^id.w" - d'w,) = 0, (2.28) 
a^J^ = 0. (2.29) 

These are in precise agreement with the relativistic formulation of Lebedev, Khalatnikov 
and Carter 0, 0, 0] provided that has the explicit form^ 



w 



i[£+p-e(/x + v- V^f] V - - [n - e(/i + V ■ Vv^)] . (2.30) 



In this approach, it is very easy as well to derive the equations of the nonrelativistic 
two-fluid model. It suffices to require the proportionality between the momentum density 

and the flow of the U{1) conserved charge, g"^ = m{nv^ + A*), where m is the mass of the 
only species of particle in the fluid. The consequences of this assumption in the context of 
effective Lagrangians for superconductors have been analyzed in Ref. 3]. By introducing 



the superfluid density = and the superfluid velocity v\= m ^dicp, one finds the same 
constitutive relations as those in Ref. |12l] : 

= p6'^ + m{n - ns)v'v^ + mnsvlvi, (2.31) 
ji = nsifi + mw ■Vs){vi-v') + {e + p)v\ (2.32) 
= ns{vl-v'). (2.33) 

Note that in this case, the particle current flow is itself conserved. 



III. SUM RULES AND HYDRODYNAMIC FLUCTUATIONS 

Here we will focus on the dynamics of fluctuations around the equilibrium state as seen 
from the rest frame of the superfluid. The velocity v(x, t) and Vy? will be non-zero only 
due to a departure from the equilibrium state. We apply the memory function formalism. 



quickly reviewed below, by following the treatment of Forster for ^He. The linearized 



^ I thank Dam T. Son for pointing out to me the overlap of these results with his unpublished work on the 
Poisson bracket approach to hydrodynamics [i, 
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hydrodynamic equations including dissipative effects will be derived as a partial check of 
the previous results. 

The fundamental quantity in this discussion is the matrix of complex response functions 

Xab{Kz)= / , (3.1) 



which is analytic for Imz ^ 0. For z in the upper half-plane, XA_B(k, + ie) is 
the Fourier transform of the equilibrium expectation value of the retarded commutator 
z6'(t)([y4(x, t), 5(0, 0)]), and the spectral function pyiB(k, i^) is the Fourier transform of 
([y4(x, t), 5(0, 0)])/2. This quantity is either real and symmetric in A ^ 5, or imaginary 
and antisymmetric. If the operators A and B have the same (opposite) signature under time 
reversal, pAB(k, tu) is odd (even) in oj. Such a signature is +1 for the energy and charge 
densities and —1 for the momentum densities and the Goldstone mode. 

When the system is perturbed by turning on a time-dependent Hamiltonian depending 
on some set of small external forces coupled to the operators {^(x, t)} 

5i/""(t) = -J2 [ t^^xA(x,t)(5F^(x,t) , (3.2) 

the standard techniques of first-order perturbation theory produce the induced changes 
5(A(x, t)). In order to follow the relaxation of the induced quantities it is convenient to 
apply external fields that are held constant for negative times and are suddenly switch off 
for positive times 

6FA{^,t) = 6FA{^)e''9i-t) , (3.3) 

where e is a positive infinitesimal number. After elimination of the external fields in favour 
of the static susceptibilities, Xab(M) = Xab(M,z = ie), through 

S{Aik,t = 0)) = J2xABik) 6FBik) , (3.4) 



the Laplace-Fourier transform of the induced changes ^(^(x, t)) becomes [14 
6{Aik,z)) = 5^-(xAB(k,z)xBc(k)-M^(^(k,i = 0)) 



BC 



(3j2CAB{k,z)xBU^)S{C{k,t = 0)) . (3.5) 



BC 



This a fundamental result that solves the initial value problem in terms of response functions. 

General arguments based on the memory function formalism make it possible to write 
the matrix C in the form 



Cab(Kz) = (3 



-1 



oo 



X(k) , (3.6) 



z-u{k)x-\k) + ia{k,z)x-\k) 
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where the matrices uj and o are given by 

u)AB^ = [ — pyii?(k,t^), (3.7) 

J-oo ^ 



oo 
oo 



aAB[k,z) = / — , (3.8) 

2m UJ ~ z 



— oo 



and they encode reactive and relaxations properties of the system, respectively. The spectral 
density 7yiB(k, cij) has the same symmetry properties as those of the matrix ujpab{^,^) 
and must define a positive quadratic form [iJ] from the requirement of positive entropy 
production. For 2; = 0, (TyiB(k, 2;) corresponds to a matrix of transport coefficients. Thus, 
in this framework, the input of ^^^(k), ^^^^(k, z = 0) and ^^^(k) for small k is all that is 
required in order to extract hydrodynamic correlation functions. With these definitions, the 
linearized hydrodynamic equations adopt the form 

z5(A(k, z)) - J2 [^AB(k) - taABiK 0)] SFeiK z) = ^<5(A(k, t = 0)) , (3.9) 

B 

where {5F(k, z)} is the set of internal forces expressing the departures of the thermodynamic 
quantities from its equilibrium values. By analogy with Eq. (13. 4p they are defined for t > 
by 

5F^(k,t) = 5^XAB(k)Wk,t)). (3.10) 

B 

To proceed further, we need some sum rules determining the appropriate x(k) and (x'(k). 
These have been derived in Ref. [13] and, essentially, they do not differ from the correspond- 
ing ones in nonrelativistic superfluids [ij] . We present the list of the required equations. 

The split of the momentum density g into irrotational gL and solenoidal gT parts leads 
to the separation of the momemtum density response function ^^^^(k, u;) = XL{k,u)k k"" + 

XT{k, uj){6^^ — k^k^) into longitudinal and transverse pieces. For a normal relativistic fiuid, the 
momentum density is given by g = /iv and the momentum susceptibility is dg^/dv^ = 6^^h 
where h is the enthalpy density. In the superfiuid phase the momentum density acquires an 
extra irrotational contribution g^ oc Vy? 

g = {h - + g,, (3.11) 

where hg is the superfiuid enthalpy density. The superfiuid momentum density arises from 
long-range order due to the Goldstone mode whose correlation function behaves as 

J~oo UJ ik^ 

where the constant ^ is positive. Thus the momentum susceptibility has a normal (isotropic) 
contribution {h—hs)S^^ and a superfiuid contribution h^k'^k^ as k —>■ 0. From the arrangement 
{h — hs)6^^ + hg^y = + {h — hs){S^^ — yy) we can write the following sum rules 



00 



limxL(fc) = lim / = (3.13) 



limxT(fc) = hm / = h-hs. (3.14) 



k^O k^O _^ TT UJ 
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Other sum rules involving if are 



^eH>{k) = [ —pe^{k,uj) = in, (3.15) 



oo 
oo 



J-oc vr 

where Eq. (13.1 5p is a consequence of the averaged Heisenberg equation of motion 

-idt{ip{^,t)) = {[H,ip{^,t)])=^^l, (3.17) 
and Xg(^(k) follows from energy conservation. The derivation of Eq. (13.151) can be made more 



rigorous by directly averaging on a restricted ?7-ensemble [l2| appropriate to superfluids'^. 
These results together with Eq. (I3.12p lead to the proportionality constant between and 
Vy?. As Xg^gs ~ Xgg, = hgk^k^ , consistency requires that 

g, = e/iV<^, h, = ^ft^ (3.18) 

in precise agreement with the linealized expression of the momentum density in Eq. (12.211) . 
Note that the quantity = fi~^'V(p plays the role of the conventional superfluid velocity. 
On the other hand, Eq. (13.131) and energy conservation yield 

r —p,l{k,u;) = hk^ (3.19) 

J —oo ^ 

omitting terms of higher order in k. Also, the equal time conmutators 

[n(x,t),g(y,t)] = tn{y,t)Vy6{^-y), (3.20) 

Kx,t),</^(y,t)] = ^5(x-y), (3.21) 

produce 

— p4(k,c.) = nfc^ (3.22) 

■oo 71" 

f°° duj , , , 

/ — pn^[k,uj) = i. (3.23) 

J-oo 

All these sum rules are valid irrespective of whether the system with broken symmetry 
is relativistic or not. In the nonrelativistic case where the mass density is not included in 
the energy density, the right side of Eq. (13.151) must be replaced by ip^R, according to the 
usual definition p = m + /xnr- Note, however, that in Eq. (I3.16P the chemical potential 
must be replaced by its leading part m. This gives the correct values appropriate to the 
nonrelativistic superfluid, hg = mris, h = mn and ^ = rig/m, where Ug is the superfluid 
particle density. 



Details of a similar computation can be fomid in p. 237 of Ref. 14 1 



8 



IV. CORRELATION FUNCTIONS AND KUBO RELATIONS 



With these results in hand, and the definition gi^(k,t) = g(k, t) ■ k, we can write the 
matrices x ^ corresponding to the four hydrodynamic variables {e, gi^, n, ip} but the 
subsequent formulae will be considerably simplified if in place of the particle density n we 
introduce a new variable q defined by 



dq = Tnd ( — 



de dn. 

n 



(4.1) 



In the linearized theory this quantity corresponds to the combination £(k, t) — {h/n)eqn{'k, t) 
representing the density of heat energy. The simplification arises because the term Ugg^ 



becomes zero, 



h/nuo^ 



0, and the only non zero ci;-terms involving q are cu, 



-isT/n as follows from the Gibbs-Duhem identity. Therefore, in terms of the 



hydrodynamic variables 



and two spectral functions 



{OA%t)} = {e,gi^,q,Lp}(k,t), 



(4.2) 

(4.3) 
(4.4) 



real and symmetric, imaginary and antisymmetric respectively, the required matrices x and 
u are given by 



/ Xe 



XAB{k) 





Xeq 







h 





Xeq 



Xqq 





\ 

T 


1 



(4.5) 



and 



lxlAB{k) 



hk 











:Ts 



—I 




ip 


Ts 



\ 



(4.6) 



/ 

hk 



where the remainder vanishing matrix elements are due to time reversal symmetry. The 
susceptibilities may be obtained from the thermodynamic potential f2 = —Vp{T,fi) by 
differentiation. These matrices are completed with the results for the transverse momentum 
density, 



x'^ 



0. 



k'k^), 



(4.7) 
(4.8) 
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From the thermodynamical derivatives 



dp 
Ye 



hXqq 



s/n 



9 ' 

eeAqq Aeq 



dp 
dq 

V dq A 



eq 



XV — V 
eeAqq Aeq 

h Xeq 

XeeXqq X 



eq 



h 



Xe 



XeeXqq Xeq 

and Eqs. fl3.1ip . (13.181) and (14. 5p we obtain the internal forces easily 

= ^ 
SF, = V, 

^(-VV + AiV- v). 



(4.9) 
(4.10) 
(4.11) 
(4.12) 



5F, 



(4.13) 
(4.14) 
(4.15) 
(4.16) 



Note that these last results depend on the choice of the linearized hydrodynamical variables. 
So, in terms of {q, gi^, n, {p} the internal forces conjugate to energy density and particle 
density are 5Fe = T~^6T and 5F„ = T6{fi/T). 

Now, we look at the memory matrix ctab (k,0) to lowest order in k. The transverse 
memory function reads 



a;^,^(k,0) = r/A;2(5*^-FA;^), 



(4.17) 



where t] > is the shear viscosity. With the previous notation, the remainder elements can 
be parametrized by 



crAB(k, 0) 



/O \ 

{C2 + lv)k^ iCik 

kTA;2 

\0 -iCik Cs / 



(4.18) 



in terms of the thermal conductivity k and three non negative longitudinal viscosities Ci, C2, 
with the constraint 

C2C3>Cl, (4.19) 

from the positivity of a. The a-matrix has been chosen with the requirement that all 
elements involving the energy density vanish at leading order in k. This is the main differ- 
ence with the corresponding matrix for the nonrelativistic superfluid in terms of variables 
{g,5(L,n,v5} 

f nTk'^ \ 

(C2 + k^ iCik 


\ -iCik Ca / 



^Ai(k,0) 



(4.20) 
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where the vanishing elements are those involving the particle density n. The reason for these 
choices lies in the fact that when a current flow is itself a conserved quantity, the dissipative 
lowest order contribution to the constitutive relation vanishes It is worthwhile to 

remark that the absence of transport coefficients for the energy current flow, Jl = g^, is only 
consistent with the usual relativistic dissipative fluid theory of Landau and Lifshitz and not 
with that of Eckart. 

All the hydrodynamic correlation functions "XabO^^^^) can be now derived from Eq. fl3.6p 
with similar results to those reported in Ref. [12] for the nonrelativistic case. They have 
poles at the lower half-plane when their denominator 

A{k, Lj) = {u^ - clk"^ + ik'^Diu){uj^ - cjk^ + ik'^D2Uj) (4.21) 

vahishes, except for the transverse correlation function with a diffusive pole at = 
—ir]k'^/{h — hs). In terms of the square of the velocity of adiabatic sound waves 

2 f dp\ hxqq 



^ =(7f = 2 ^ (4-22) 

and the normal enthalpy density hji = h — C.fi'^, the speeds of propagation of flrst and second 
sound are given by 

cl + cl = c^fl + ^^^V (4.23) 



hnXqq 

clcl = (4-24) 

n^KXqq 

and the attenuation constants can be written in the form 

1 4 \ f(/lC3-2uCl) C^TXeeK 

hn \ 3 y /l„ hXqq 

2n , 2n _ c'an'fi\qq + 2nsTfiXeq + s'T\ee) , 4 



1T?hKXqq V 3 

c^e/iCs 2c2e(n/ix,g + sTxeq)Ci 



/ijj Tlh^Xqq 

A (4.26) 

Xqq 



Some Ward identities can be directly checked. We list them 

^Xee{^.^) = kxg^eik,uj), (4.27) 

^XegAk,^) = kxLik,uj) - hk, (4.28) 

UJXenik,Uj) = kxg^nik,uj), (4.29) 

uJXev>ik,uj) = kxg^^ik,uj)-iij,. (4.30) 

However, Ward identities for correlation functions involving currents such as t^' or j* cannot 
be checked within this approach because these flows are not hydrodynamical modes. 
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lim lim — Im Xee{k, to) = lim lim — Im Xeqik, cu) = 0. (4.36) 



The usual catalog of Kubo relations giving the transport coefficients can be obtained 
from the following limits 

f] = \im\im — lmxTik,uj), (4-31) 

4 uj 

+ = linilim — ImxL(fc,^), (4.32) 
3 t^^o k^o 

'^^ r = limlim — Imx„n(^,tu), (4.33) 

Ci = - lini lim^ReXgL¥'(^'^)' (4-34) 
Cs = lim lim culm x<z, a, (A;, t^), (4.35) 

where for the thermal conductivity we have replaced Xqq by {h/riYxnn since 

7^ Im Xeeik, u) = lim lim ^ 

The shear and bulk viscosities 77 and (2 of the superfluid quark matter in the color-flavor 
locked phase have been computed in Refs. 15|, [l6| , and all the bulk viscosities in neutron 
stars have been computed in Ref. [l7| . 

Finally, the constitutive relations from Eq. (13.91) read 

Hie) = (/i-e/i')v + e/iVy., (4.37) 

S{in) = in- e/.)v + ^V^-K V (I) , (4.38) 

V*''(x,t) = 5/i(x,t) - CiV • V - ^3 (VV - /^V ■ v), (4.39) 

5{f^) = 5p{x, t)5'' - ri (^ViVj + V^v^ - ^ V • v 6'^ 

-5'^ (C2V ■ V + eCi(VV - AiV ■ v)) . (4.40) 

The reactive parts of these equations agree to linear order with those found in Sec. HTl 

Summarizing, we have presented the distinctive features of relativistic superfiuids (in 
comparison with those in the nonrelativistic regime) accounting for their hydrodynamics and 
the correlation functions have been derived. The crucial assumption is that the momentum 
density coincides with the flow of the conserved energy. As a consequence, there are no 
dissipative contribution to the energy current and the thermal conductivity appears in the 
current of the conserved particle number, according to the fluid theory of Landau and 
Lifshitz. 
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